The chiral representation of the -kN scattering amplitude and the pion-nucleon sigma term 
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We present a novel analysis of the -kN scattering amplitude in Lorentz covariant baryon chiral perturbation 
theory renormalized in the extended-on-mass-shell scheme. This amplitude, valid up to 0{p'^) in the chiral 
expansion, systematically includes the effects of the A(1232) in the (5-counting, has the right analytic properties 
and is renormalization-scale independent. This approach overcomes the limitations that previous chiral analyses 
of the TiN scattering amplitude had, providing an accurate description of the partial wave phase shifts of the 
Karlsruhe-Helsinki and George-Washington groups up to energies just below the resonance region. We also 
study the solution of the Matsinos group which focuses on the parameterization of the data at low energies. 
Once the values of the low-energy constants are determined by adjusting the center-of-mass energy dependence 
of the amplitude to the scattering data, we obtain predictions on different observables. In particular, we extract an 
accurate value for the pion-nucleon sigma term, ct^jv. This allows us to avoid the usual method of extrapolation 
to the unphysical region of the amplitude. Our study indicates that the inclusion of modem meson-factory and 
pionic-atom data favors relatively large values of the sigma term. We report the value ctttn = 59(7) MeV and 
comment on implications that this result may have. 

PACS numbers: 13.75.Gx, 11.30.Rd, 12.39.Fe. 13.85.Dz 



The sigma terms, cTt^n and CTs, are observables of funda- 
mental importance that embody the internal scalar structure 
of the nucleon, becoming an essential piece to understand the 
origin of the mass of the ordinary matter. The pion-nucleon 
sigma term, (Jt^n, is a key ingredient in investigations of the 
QCD phase diagram and in the study of nuclear systems 
On the other hand, cJt^m and (Js appear as hadronic matrix el- 
ements in the neutralino-nucleon elastic scattering cross sec- 
tion. Unfortunately, our current knowledge of the sigma terms 
is far from satisfactory. With the advent of experimental re- 
sults on dark-matter searches, different authors have pled for 
a more accurate experimental determination of these quanti- 
ties 13-5]. 

The CT^AT is defined as the nucleon matrix element of the 
light-quark scalar current. 
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where rh — (to„ + fnd)/2, with m„, rrid the light-quark 
masses. The sigma term can be obtained from the ttA^ scat- 
tering data by extrapolating the scattering amplitude to the 
Cheng-Dashen point which lies in the unphysical re- 

gion of the Mandelstam plane. The usual method to perform 
this extrapolation is by means of an energy-dependent param- 
eterization of the data in partial waves (PW) supplemented by 
dispersion relations that impose strong analyticity and unitar- 
ity constraints onto the scattering amplitude at low energies. 
The current uncertainty in a^^N originates from discrepancies 
between the classical PW analysis of the Karlsruhe-Helsinki 
(KH) [18] group and the more modern one performed by the 
George- Washington |9] (GW) group. More precisely, the KH 
amplitudes were used b y G asser et al. to obtain the canonical 
result (TtAt — 45 MeV lUOll . whereas the analysis of the GW 
group, which includes modern meson factory data, leads to a 
larger value ct^at = 64(7) MeV 10]. 

The main difficulty of the traditional method to obtain (TttW 
is assessing the errors that propagate in the extrapolation to the 



Cheng-Dashen point from the systematic uncertainties associ- 
ated to a particular parametrization of the data. These prob- 
lems, together with the persisting discrepancy in the values 
reported by the different PW analyses, have led to new strate- 
gies for the determination of the sigma terms. Particularly 
noteworthy is the intense campaign developed by the LQCD 
community to calculate these matrix elements using new pow- 
erful algorithms and computational resources llT2ll . 

In this paper, we focus on the extraction of the a-j^N 
from TriV scattering data and using chiral perturbation theory 
(xPT), which is the effective field theory of QCD at low en- 
ergies ll^l - (T6ll . This is a suitable framework to shed light on 
the experimental discrepancies since it allows for an investi- 
gation of the chiral Ward identity that relates the isoscalar ttN 
scattering amplitude and cr^Ar, giving a handle on the eiTors 
committed at each order of the power counting. In fact, one 
recovers the low-energy theorem at the Cheng-Dashen point 
that is exploited by the dispersive methods mentioned above. 

In xPT, one can alternatively use a more elegant manifes- 
tation of the same Ward identity between the two observ- 
ables fT?]. At 0{p^), cr^jv only depends on one apriori un- 
known low-energy constant (LEC), ci (cf. Eq. (2) below). 
Due to the non-linear realization of chiral symmetry under- 
pinning xPT, this LEC also contributes to nucleon processes 
with an even number of external pion legs and, in particular, 
to the isoscalar part of the TriV scattering amplitude. There- 
fore, determining the value of this constant with a fit to the 
scattering data allows to predict (TttW, avoiding any analyti- 
cal extrapolation of the amplitude onto the unphysical region 
(as the one to the Cheng-Dashen point used in the dispersive 
analyses). 

The low-energy structure of the ttN scattering amplitude 
has been studied within different approaches tackling the sub- 
tleties in the power counting that appear in the baryon sec- 
tor of xPT (BxPT) (for reviews see Refs. BHltl). After 
the seminal paper of Gasser et al. ITtIi . it was first studied in 



TABLE I: Physical observables obtained from the 0{p'^) tvN scat- 
tering amplitude in the EOMS renormalization scheme fitted to dif- 
ferent PW analyses up to Wmax = 1.2 GeV (W„ax = 1.16 GeV 
for EM). The scattering lengths are in units of 10^^ wi^^. 





Xd.o.f. 


hA 


9ttN 


Agt [%] 


«o+ 


«o+ 




0.75 


3.02(4) 


13.51(10) 


4.9(8) 


-1.2(8) 


8.7(2) 


GW|2] 


0.23 


2.87(4) 


13.15(10) 


2.1(8) 


-0.4(7) 


8.2(2) 


EM [39] 


0.11 


2.99(2) 


13.12(5) 


1.9(4) 


0.2(3) 


7.7(1) 



heavy-baryon (HB) xPT fl?) by Fettes et al. up to 0{p^)\\^ 
and 0{p^) |20] in the chiral counting. In these works, a pre- 
cise description of the PWs was obtained at low energies, al- 
though the values of the LECs contain important contributions 
from the A(1232) resonance and the results for the aT^M were 
not accurate, being typically too large. The inclusion of the 
A as an explicit degree of freedom in the so-called small- 
scale-expansion (SSE) ll2Tll (that counts e = A/a — Mjy ^ p) 
up to 0{e^) 1 22], offers a noticeably increase in the range 
of energies described compared with HBxPT at 0{p^) ||19|]. 
Nonetheless, there is a strong dependence on the fitted values 
of the LECs with the PW analysis used as input that prevents a 
direct extraction of aj^N by fitting scattering data |22]. After 
these difficulties, the conclusion was that the chiral conver- 
gence was not fast enough in the physical region so to extract 
useful information on ct^^t from the PW phase shifts [23]. 

It has been shown that the non-relativistic expansion im- 
plemented in the HB approach does not converge in part of 
the low energy region Il4l24l [25ll. This led to the studies in 
the manifestly Lorentz covariant infrared (IR) BxPT Ii24 l26l - 
I28I1 . In this case, the amplitude up to 0{p'^) without the A 
as explicit degree of freedom, shows an accurate and rapidly 
convergent description in the subthreshold region but fails to 
connect it to the physical one [26], confirming the conclusions 
about CTjrJV drawn from the previous works in HB. 

In this paper, we present a xPT analysis of the ttN- 
scattering amplitude and of the pion-nucleon sigma term 
up to 0{p^) accuracy that includes two main improve- 
ments over previous work. In the first place, we use 
Lorentz covariant BxPT with a consistent power counting ob- 
tained via the extended-on-mass-shell (EOMS) renormaliza- 
tion scheme |29]. This prescription, instead of IR, is used be- 
cause the latter introduces unphysical cuts that may influence 
the low-energy region [ 16l l30l |31|] . As it has been recently 
shown in Ref. [32], the huge Goldberger-Treiman (GT) dis- 
crepancy, of ^ 20%, found in the IR scheme [|27i I28II can 
be traced back to the analyticity issues of this method rather 
than to a breaking of the chiral convergence in the ttN sys- 
tem. In addition, we obtain amplitudes independent of the 
renormalization scale, which is not the case for those given by 
IR [24, 26-28]. Secondly, we explicitly include the A taking 
into account that, below the resonance region, the diagrams 
with the A are suppressed in comparison with those with the 
nucleon. This can be implemented in the so-called (5-counting 
by assigning an extra fractional suppression of 0{p^^^) to the 
A-propagators in the Feynman diagrams ll33ll34ll . For the A^A 




FIG. 1: (Color on-line) Phase shifts given by the Lorentz covariant 
O(p^) nN scattering amplitude in the EOMS scheme fitted to the 
GW solution (circles) |9] up to W,nax = 1.2 GeV. 



chiral L agrangi ans, we use the consistent formulation of Pas- 
calutsa Il34l - [36ll which filters the unphysical components of 
the relativistic spin-3/2 spinors and eliminate the dependence 
on off-shell parameters that the conventional vertices have. 
The technical details of this calculation and the complete re- 
sults derived thereafter are presented in detail elsewhere ll32ll . 
In the following, we outline the analysis and show its main 
results on the vrA^ phase shifts and, more specifically, on the 
pion-nucleon sigma term. 

The calculation proceeds as in Ref. ll28ll but with the loops 
of that reference treated in the EOMS scheme. This is 
achieved by canceling the UV divergences obtained in dimen- 
sional regularization such that the power-counting breaking 
pieces of the loops are absorbed into the 0{p) LECs, gA (axial 
coupling of the nucleon) and M^, and into the 4 0{p^) LECs, 
ci, C2, C3, C4. The 5 combinations of 0{p^) LECs, di + d2, 
da, c?5, c?i4 — (ii5 and dig are renormalized in the MS scheme. 
Besides that, we also include the Born-term with an interme- 
diate A(1232) resonance and leading 0{p) vertices given by 
the A axial coupling Ha- The Born-terms with 0{p^) iVA 
couplings ll22l[37ll38ll have also been considered but they give 
a negligible contribution and have been omitted in the present 
study, whereas the coiTesponding loops with A propagators 
are of higher-order. 

We fix the values of the LECs fitting the center-of-mass 
(CM) energy dependence of the 2 S- and 4 P-wave phase 
shifts obtained from the chiral amplitude to the latest solutions 
of the KH [8] and GW |9] groups. In addition, we include 
the analysis of the Matsinos' group (EM) [39] which focuses 
on the PW parameterization of the data at very low energies 
without imposing dispersive constraints from the high-energy 
region. We follow the logic of Ref. |28] to assign errors to the 
first two analyses (they do not provide eiTors) while for the lat- 
ter we include the eiTors provided there. The fits are done from 
the lowest CM energies above threshold, Wth — 1.08 GeV, 
up to Wmax = 1.2 GeV which is below the A(1232) re- 
gion (the EM analysis only reaches Wmax — 116 GeV). 
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The parameters cja, Mjv, Ma, ra-^ and f.^ are fixed to their 
experimental values |28|]. Although the A^A axial coupling 
can be determined using the A(1232) width, we also fit this 
constant to the PW phase shifts. The suitable value to com- 
pare with is the one obtained from the Breit- Wiener width 
Fa = 118(2) MeV [40], namely Ha = 2.90(2) H. 



TABLE 11: Values of the Oip^) LECs in units of GeV"! and of g^m 
in MeV obtained from the different -kN PW analyses. 





Cl 


C2 


C3 


C4 


O'wN 


KH 


-0.80(6) 


1.12(13) 


-2.96(15) 


2.00(7) 


43(5) 


GW 


-1.00(4) 


1.01(4) 


-3.04(2) 


2.02(1) 


59(4) 


EM 


-1.00(1) 


0.58(3) 


-2.51(4) 


1.77(2) 


59(2) 



The results on a selection of physical observables are shown 
in TableU The errors quoted there are only of statistical origin 
and additional theoretical uncertainties are to be added. In 
Fig. [U we plot the phase shifts of the S- and P-waves given 
by the vrA^ scattering amplitude in the EOMS scheme and at 
0{p^) in the (5-counting, fitted to the GW solution (circles) [9] 
up to Wmax = 1.2 GeV. Similar plots can be obtained for the 
KH and EM solutions. The figure shows that the description 
of the lowest PWs is very accurate up to energies below the 
A-resonance region, covering a range of energies larger than 
in previous perturbative analyses |l 20il22ll27il28ll . The quality 
of the description is reflected by the small q f listed in the 
second column of Table |I] which furthermore shows that the 
description of GW and EM PW analyses is better than the KH 
one. 

As one can infer from the third column of Table |I] only the 
GW solution gives a result on Ha that is perfectly compatible 
with the determination from the A width. In the fourth col- 
umn, we show the values obtained for the TriV coupling that, 
compared with the axial coupling gA, gives the GT discrep- 
ancy Act in the fifth column. These can be compared with 
the numbers independently extracted from A^A^- scattering 
{g^N ^ 13.0) [41] and pionic atom {g^N = 13.12(9)) ll42n 
data. For the discussion of the large violation obtained in 
IR ByPT IztI I28I1 as compared with the one obtained here, 
see Ref. Il32ll . Results for the isoscalar (a^^) and isovector 
(cq^) scattering lengths are shown in the last two columns 



of TableU Minding that changes of 5-10% can be easily ex- 
pected from higher-order and isospin corrections [430, we 
can compare with the values independently extracted from 
pionic-atom data, a^^ = 0.0861(9) and a'^^ 0.0076(31) 
TO^^ [42]. The impact that the pionic-atom result for aQj_ has 
on the value of a-^^ lfloll44ll is addressed below. 

As mentioned above, the isoscalar TriV scattering amplitude 
is related with a-^^ through the LEC ci . In the first columns of 
Table HH we show the fitted values for the C(p^) LECs. The 
values and errors quoted there correspond to the mean and 
standard deviation obtained after considering fits to the KH 
and GW PW phase shifts for various Wmax, from 1.14 GeV 
to 1.2 GeV {Wmax ^ 1-16 GeV for the EM analysis) in in- 
tervals of 0.01 GeV. The purpose of this strategy is to take 
into account the dispersion of the values of these LECs (and 
of (Jtjm) against the data set included in the fits. As we can 
see, our results remain stable to the increase of the maximum 
energy and to the particular analysis used as experimental in- 
put. This is in remarkable contrast with the strong sensitivity 
of the values of the LECs obtained in the 0{€^) study done in 
HBxPT-SSE ii. 

On the other hand, the results in Table are quite dif- 
ferent to the ones obtained without the explicit inclusion of 
the A. In this case, we can describe the GW phase shifts 
up to Wmax ^ 1.14 GeV ixlof = 0.62) and we ob- 
tain Cl = -1.54(5), C2 = 3.92(6), C3 = -6.87(6) and 
C4 — 3.79(3) (all in units of GeV~^). Comparing with the 
values in Table [III we see that the contribution of the A to 
the 0{p^) LECs C2-4 is compatible with the one estimated 
by resonance saturation hypothesis ||45|] . = 1.9... 3. 8, 
= -3.8 ... - 3 and = 1.4 .. . 2.0 (in GeV"^ For 
the Cl counter-term the A contribution is neglig ible 1^, 45!|. 
We interpret the difference, of around 0.5 GeV^^, between 
our result in the second row of Table |ll] including the A and 
that without this resonance, as a clear indication that the LECs 
are stabilized once the tree-level A exchange contributions are 
taken into account llssi l46tl . 

We calculate Gt^n at 0{p^) employing covariant BxPT in 
the EOMS renormalization scheme. The pion-nucleon sigma 
term can be obtained either from the scalar form factor of the 
nucleon, Eq. ([T]i, or from the quark mass dependence of its 
mass and the Hellmann-Feynman theorem. The resulting ex- 
pression is 



cj.N = -4ci774 - I ^ \ arccos — ^ + log — ^ , (2) 



and leads, in the non-relativistic limit, to the HB result up to 
0{p^) (aT^ . Besides the error propagated from ci, this ex- 
pression carries a theoretical uncertainty coming from higher- 
order contributions. We estimate this by computing the next 
subleading correction, at 0(p^/^) in the J-counting, which is 



given by a loop diagram with an insertion of a A propaga- 
tor [48]. This amounts to a contribution of —6 MeV (that we 
take as an irreducible uncertainty in our determination) to be 
compared with the one at 0{p^) of —19 MeV. Furthermore, 
we have calculated the 0{p'^) corrections given by the loop 
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diagrams with an insertion of the 0{p^ ) LECs f^T). With the 
values of the LECs in Table HIl we obtain that they span from 
—2 to —4 MeV. These results suggest a clear convergence pat- 
tern for the chiral expansion of (Tttn as well as they confirm 
the hierarchy at low energies between the nucleon and A con- 
tributions that is implemented in the (5-counting. 

With Eq. (|2]l and the values for ci obtained in the inter- 
val 1.14 GeV < Wmax < 1-2 GeV, we determine the means 
and standard deviations of a^rN listed in the last column of 
Table Hn The values of ct^^t extracted from the different anal- 
ysis are not completely consistent among each other The KH 
number reproduces the canonical result cr^jv — 45 MeV lHoll . 
whereas those determined from the GW and EM solutions 
agree with the dispersive result of the GW group cr^jv = 
64(7) MeV 1 1 1]. Furthermore, the result from the EM anal- 
ysis also agrees with a-r^N — 56(9) MeV, obtained by 01s- 
son iiill using a dispersive sum-rule and the threshold param- 
eters provided by an early version of the EM solution. 

Although our results for each of the PW solutions are 
consistent with those obtained extrapolating the data to the 
Cheng-Dashen point, the B^PT approach applied here relies 
solely on the information in the region where the data actually 
exist. We then give an estimation on the uncertainty commit- 
ted in the relation between (Jt^n and the ttN scattering ampli- 
tude that is based on effective field theory grounds. On top 
of that, the dispersion of the results with respect Wmax and 
among the different analyses allows to disentangle the sys- 
tematics coming either from the data basis employed or the 
particular parameterization of the data. In this sense, the con- 
sistency between the results derived from the GW and EM so- 
lutions is very remarkable since these are quite different PW 
parameterizations having both in common the inclusion of the 
wealth of low-energy data collected along the last 20 years 
in meson and pion factories lUi [IK] with many points not in- 
cluded in KH 1 8]. Therefore, our results suggest that the mod- 
ern meson scattering data lead to a value for cr^^y larger than 
the one obtained from the older KH analysis IllOtl . Neverthe- 
less, a re-analysis of the modern data set with the KH method 
would be extremely valuable in order to reach a definite con- 
clusion in this regard (see e.g. Ref. |49]). 

Another important and independent source of information 
comes from the pionic-atom data on Qq^. It has been noted 
before in dispersive studies I^Tol [Tlll44ll that the sign of this ob- 
servable is strongly correlated with the value of CTjt at ■ While 
the KH result is compatible with the old negative results, it 
is not anymore with the recent positive values extracted from 
modern pionic-atom data and using improved phenomenolog- 
ical approaches |42]. These are, on the other hand, compatible 
with the scattering data determinations obtained from the GW 
and EM solutions. The effect that a non-negative result on 
aj_i_ has on a-^N was quantitatively studied by the GW group 
concluding that a value of a^J^ > produces a raise on the 
sigma term of, at least, 7 MeV [11]. 

Finally, we want to emphasize that only our results based 
in the GW analysis are perfectly compatible with all the phe- 
nomenology that can be extracted from independent experi- 
mental sources. We remind here that the KH analysis gives 
rise to a value for Ha that is not compatible with the value 



obtained from the A(1232) width (in agreement with the KH 
overestimation of this observable) and to a value for g^^N that 
leads to a sizable violation of the GT relation, which is nowa- 
days theoretically implausible. As for our study of the EM 
PW analysis, we found a value for the isovector scattering 
length that is too small as compared with the accurate values 
obtained from pion-atoms data f3^j4^. 

With these considerations, one obtains the following value 
for aTrN, as it is extracted from the analysis of ttN modern 
scattering data 1^, ^] and using Lorentz covariant B^PT in 
the EOMS scheme up to 0{p'^) in the (5-counting, 

cr^^ = 59(7)MeV. (3) 

The error includes the higher-order uncertainty estimated 
above added in quadrature with the one given by the disper- 
sion of the values in the average of the GW and EM results. If 
one were to include the KH result in this estimation, the result 
would be slightly reduced by 2-3 MeV. 

As a concluding observation we want to address the fact 
that this relatively large value of (TttW may appear to be in 
conflict with some established phenomenology. In particu- 
lar, it may give a new twist to the old puzzle concerning the 
strangeness content of the nucleon lUlll . This is based on 
the relation that is obtained in HB^PT up to ©(p^) accuracy 
among the SU (3) i?-breaking of the baryon-octet masses, aT^N 
and the observable y quantifying the strangeness content of 
the nucleon ifsollsill . For the value of the sigma term obtained 
in the present work, this relation leads to a contribution of the 
strange quark to the nucleon mass of several hundreds of MeV. 
It is interesting to note that the usual method to derive this re- 
lation does not include explicitly the effects of the decuplet 
resonances, which have been shown to largely cancel those of 
the octet in the strangeness content of the nucleon ll47i l52ll . 
This is, indeed, consistent with recent B^PT determinations 
of the sigma terms using LQCD results on the baryon masses 
and explicitly including the decuplet contributions jssl Isill . 
showing that a relatively large value of a^^N ~ 60 MeV ll54ll 
is not at odds with a negligible strangeness in the nucleon. 
Another caveat arises in chiral approaches to nuclear matter, 
in which a large value of o-ttTV would lead to a vanishing quark 
condensate at too low densities It is important to note 

that a non-zero value of the in-medium temporal component 
of the pion axial coupling, ft, is also a necessary condition for 
the spontaneous breaking of chiral symmetry ]55]. Hence, an 
analysis of the density dependence of this quantity, together 
with the quark condensate, is necessary in order to properly 
discuss about chiral symmetry restoration in nuclear matter 

In summary, we have presented a novel analysis of the ttN 
scattering amplitude in Lorentz covariant B^PT within the 
EOMS scheme up to 0{p^) and including the effects of the 
A(1232) explicitly in the 5-counting. This covariant approach 
ensures the right analytic properties of the tree-level and loop 
corrections to the amplitude, providing a model-independent 
framework to comprehensively and accurately study the phe- 
nomenology associated with the different PW parameteriza- 
tions of the TriV-scattering data. In particular, we found that 
we perfectly describe the PW phase shifts of the KH, GW 
and EM groups up to energies below the A-resonance region. 
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at the same time as we agree in the values of the scattering 
observables. It is worth stressing that, apart from the phase 
shifts, our results using the GW analysis are perfectly com- 
patible on important observables with those obtained from in- 
dependent phenomenological sources. 

We show that our amplitudes are suitable to extract an ac- 
curate value of a-i^N from scattering data and avoiding the ex- 
trapolation to the unphysical region using a method based on 
EOMS-BxPT. Namely, the pion-nucleon sigma term can be 
properly calculated at 0{p^) only when B^PT is formulated 
in a way that keeps the analytical properties of the amplitude, 
accounts for the important effects of the A resonance in the 
LECs and gives results independent of the renormalization 
scale. It follows that the extraction method to calculate a-j^N 
is under good theoretical control. Consequently, we ratify the 
discrepancy between the KH and GW groups and give support 
to the latter, which is in agreement with the one that we obtain 



from the study of the latest EM solution. We conclude that 
recent analyses of the modern data lead to a relatively high 
value of (Jtjn, cf. Eq. 
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